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ABSTRACT 

An assay technique for quantifying 
milk-clotting enzymes is described. By 
the technique, enzyme concentrations as 
small as 2 x 10 — 5 rennet units/ml can be 
determined accurately in milk solutions 
in as little as 2 h. The technique can be 
used to assay for enzymes in milk and 
whey solutions without prior removal of 
casein. The technique should find ap- 
plication in detecting enzymes in a 
variety of systems commonly encountered 
in studies of milk and cheese products. 

INTRODUCTION 

Quantification of milk-clotting activity in 
solutions containing proteolytic enzymes is a 
major concern in industrial cheese making 
and cheese research. Many assay techniques for 
quantifying renneting activity have been 
described (1, 3, 4, 5, 9, 11, 14, 15, 19, 20, 
22, 23, 24, 25, 28). Immunological agglutination 
reactions have been used to detect chymosin, 
porcine pepsin, and bovine pepsin in cheese and 
milk products and to make quantitative deter- 
minations of the enzyme content in rennet 
solutions (9, 24). Several assay techniques 
requiring quantitative analysis for a product of 
the hydrolysis reaction catalyzed by the enzyme 
to be assayed also have been reported. Sub- 
strates that have been used to assay various 
milk-clotting enzymes are: 1) milk (4), 2) 
casein (22), 3) hemoglobin (1), 4) synthesized 
peptides (23, 25), and 5) nonpeptide substrates 
such as phenyl sulfite (28). Special preparative 
techniques, which utilize milk and casein as 
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substrates and which allow the reaction product 
to be detected using absorbance measurement 
at 217 nm (11), Kjeldahl nitrogen assays (5), 
or other protein determination techniques [e.g., 
Lowry et al., (20)] , have been used to quantify 
the product released. 

Casein-agar diffusion techniques (3, 14, 15, 
19) can be used to determine independently 
milk-clotting and general proteolytic activity 
without requiring direct analysis for a product 
of the hydrolysis reaction. These methods are 
relatively accurate and sensitive but incon- 
venient; the casein-agar must be specially han- 
dled in controlled temperature areas-, each 
preparation of agar must be calibrated with 
known enzyme solutions; and assays take up to 
48 h to complete. 

Historically, the most common method for 
determining enzymatic activity of proteases has 
involved measuring the time a milk sample 
takes to clot after addition of enzyme and 
relating this time to the quantity of enzyme 
added. The clotting point can be determined 
visually when flocculation of milk precipitate 
begins (2). It can be taken as the time when a 
rapid increase of viscosity of the milk solution 
occurs (17, 18) or as the time the milk forms a 
gel (30). Experimental data indicate that there 
is a linear relationship between clotting time as 
defined by these methods and the inverse of the 
enzyme concentration (16, 17). 

Clotting determinations using raw skim milk 
are not particularly sensitive or accurate as the 
pH of skim milk (ca. 6.7) is far from the 
optimum for enzymatic action of rennet and 
other acid proteases used for clotting milk (13). 
Even at lower pH, skim milk takes a relatively 
long time to clot. Furthermore, reproducibility 
of determinations with skim milk is questionable 
because "clottability" depends to a degree on 
the origin and history of the milk sample (33). 
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An improved substrate for assay purposes is 
reconstituted nonfat dry milk (2). Calcium 
often is added to solutions of this material 
to promote rapid clotting. Further enhancement 
of clotting rate and stabilization of the enzyme 
being assayed is accomplished by lowering the 
pH of the milk solution to around pH 6.0. 
Diluting the reconstituted milk to .1 to 5.0% 
solids provides a solution with a more distinct 
endpoint for visible coagulation than fully 
reconstituted milk. The coagulation time in 
these solutions is substantially less than in skim 
milk for a given enzyme concentration (12). 
Very dilute milk solutions (less than .5% solids) 
can be used for spectrophotometric deter- 
mination of the coagulation time (21). 

The simple "clotting time" assays described 
are limited as to the minimum enzymatic 
activity that can be detected accurately. This 
limit is caused by three factors: 

1) Long coagulation times characteristic of 
low enzyme in the assay solution result in 
bacterial degradation of milk proteins and 
changes of clotting properties during assay. 

2) Small enzyme concentrations cause the 
physical flocculation process to proceed so 
slowly that it is difficult to specify the exact 
point of coagulation. 

3) At low enzyme concentrations the rela- 
tionship between clotting time and the inverse 
of enzyme concentration is nonlinear (10). 
Hence, it is difficult to develop a simple func- 
tional equation relating these variables. 

Our paper describes a method of increasing 
sensitivity of the traditional milk-clotting assay. 
It avoids the aforementioned problems without 
greatly increasing the inconvenience of the 
technique. The method involves incubating the 
test sample in diluted reconstituted skim milk 
and adding an additional amount of enzyme to 
clot the sample. Observed clotting time can be 
related simply to initial amount of enzyme in 
the sample. 

Mathematical Basis of Assay Technique 

Flocculation of casein proteins in the 
presence of calcium ions is highly dependent on 
the relative amounts of a s - and K -casein in the 
solution (7, 32). The principle action of milk- 
clotting enzymes is to hydrolyze K -casein (31), 
thereby reducing the ratio of K :0i s in the 
caseinate system. When this ratio has been 
reduced sufficiently, a precipitate of calcium 



caseinate is seen readily if the calcium ion 
concentration in the solution is high enough. 
The time at which large agglomerates of floc- 
culated calcium caseinate become visible to the 
naked eye commonly is called "clotting time". 

The enzymatic hydrolysis of K-casein appears 
to follow Michaelis-Menten kinetics (11) i.e. 

rate of K-casein hydrolysis 

= V m (K)/[K m + 0O] J [1] 

where (k) refers to concentration of K-casein, 
V m is the maximum reaction rate, and K m is 
the Michaelis constant. 

The reaction rate also may be written: 

dK/dt = [-k cat E (/c)/[K m + (K)] [2] 

where t is the reaction time and V m has been 
replaced by the product of the catalytic con- 
stant, k cat , and Eo- The E is the total con- 
centration of enzyme (both bound and free). 

For present purposes it is sufficient that the 
rate of hydrolysis for a fixed substrate con- 
centration be proportional to the total con- 
centration of enzyme (Eo). Hence, equation [2] 
can be written in the form : 

d(K)/dt = E /f(K) [3] 

where for the case described: 

f(K) = -[K m + (K)]/k cat (K) [4] 

More generally f(K) can take on a number of 
forms characteristic of steady-state kinetics (6, 
26). 

If one assumes that no flocculation of casein 
occurs until a critical amount of K-casein has 
been hydrolyzed (7) and that thereafter clotting 
occurs immediately, the relationship between 
the clotting time (t c ) and total enzyme con- 
centration can be found by separating variables 
and integrating equation [3] from t=0; (k) = 
(K ) to t=t c; (K) = (K C ). 

(K c ) t c 

] f(K)d(K)= j E dt [5] 
(K ) 

where (K ) is the initial concentration of 
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K-casein and (k c ) the K-casein concentration at 
the moment clotting takes place. 

For a constant concentration of enzyme this 
gives: 



the secondary addition of enzyme is less than 
the clotting time when the two amounts 
are added simultaneously. In this case the 
mathematics are: 



where : 



t c = S2 C /E 
(K c ) 



I2 C = ) f(K)d(K) 

(Ko) 



16] 



17] 



Equation [6] predicts that the clotting time is 
inversely proportional to the total enzyme 
concentration. This form has been noted 
by several investigators (17, 29). 

The experimentally observed clotting time, 
t c (obs), usually is related to the enzyme con- 
centration by an equation of the following 
form (10, 16): 



t c (obs) = f2 c /E + t 



[8] 



where t is an experimentally determined 
constant. Verification of the form of equation 
[8] is shown for a limited range of enzyme 
concentrations in Figure 1. 

This correlation can be explained by assuming 
that a short time (t ) must elapse before 
clotting is observed even after the critical 
amount of K-casein has been hydrolyzed, i.e., 



t c (obs) = t c + t 



19] 



In principle, equation [8] can be used as a 
basis for measuring the clotting activity of an 
unknown sample. This is accomplished by 
measuring the observed clotting time for such a 
sample and using the parameters t and n c 
obtained from solutions of known enzyme 
concentrations to calculate the concentration 
of the unknown (17). However, this simple 
procedure becomes impractical for very low 
enzyme concentrations for the reasons men- 
tioned. The discussion that follows outlines the 
basis of a method for extending the limit of 
detection of the standard milk-clotting assay. 

If a small amount of enzyme is added to a 
casein solution first, and if at some time before 
coagulation another amount of enzyme is 
added to the milk, the clotting time following 



From t=0 to t=tj , 

f(K)d(K) = ) E x dt [10] 

(K ) 

where (Ki) is the concentration of K-casein 
remaining at the time (ti) when the secondary 
addition of enzyme is made, and Ei is the total 
enzyme concentration (bound plus free) due to 
the primary addition of enzyme. 



From t=ti to t=t c 



f(K)d(K)= J (Ei +E 2 )dt 
(Kl) tj 



[11] 



where F,2 is the concentration of enzyme 
(bound plus free) due to the secondary addition 
of enzyme. 

Addition of equations [10] and [11] gives.- 

(K c ) tj t c 

f(K)d(«)= j Eidt! + j (Ei +E 2 )dt 
(«o) t, [12 ] 



2 100 




s \= 1,6 (sec) 



25 50 75 100 

(CONCENTRATION OF ENZYME)" 1 (RU/ralT 1 

Figure 1. Experimental verification of equation 
[8] . (Assay solution clotted with rennet extract). 
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For constant Ei and E 2 integration yields: 
fi c = E 1 t 1 + (E t +E 2 )(t c -t!) 

[13] 

Again if t c = t c (obs) — t : 

D, c = Ext! + (Ei + E 2 ) (t c (obs) - ti - t ). 

[14] 

The term (t c (obs) — ti) is the observed clotting 
time following the secondary enzyme addition. 
This difference will be termed t a . Hence: 

S2 c = Eiti +(Ei +E 2 )(t a -t ) [15] 

Division by (Ei + E 2 ) and rearrangement gives: 



E l 

t = tj + + t 

Ej + E 2 (Ei + E 2 ) 

[16] 

Further, from equation [8] the final two terms 
in equation [16] represent clotting time (t c ') 
characteristic of an enzyme concentration (Ei 
+ E 2 ) added to the assay solution at time ti = 
0. Hence, 

^c 

t c = + t [17] 

(Ei+E 2 ) 

With this substitution equation [16] then 
becomes: 

t a = -- — t, + t c [18] 

Ei + E 2 

This equation predicts that the observed 
secondary clotting time (t a ) will decrease in a 
linear fashion as the time when the secondary 
enzyme is added (ti ) is increased. The slope of 
a plot of t a vs. ti is related to the relative 
enzyme concentrations Ei and E 2 by the 
equation: 

Ei r 

slope = = [19] 

Ei+E 2 (r+1) 

where r = Ei/E 2 . 



The intercept t c ' is related to the clotting time 
observed when a concentration of enzyme E 2 is 
used to clot the casein solution (t c ) by the 
equation, 

t c ' = (tc-t ) +t [20] 

E 1+ E 2 

Figure 2 shows how these variables are related. 

The relationship between t a and t! can be 
used to determine the concentration Ei in a 
casein solution. By determining the slope 
of a plot of t a vs. ti, the ratio of the two 
enzyme concentrations Ei and E 2 can be 
estimated because, 

Ei/E 2 =-slope/(l + slope) [21] 

If concentration E 2 is known, Ei can be found. 

MATERIALS AND METHODS 

Pasteurized, low temperature, nonfat dry 
milk from Consolidated Badger Co-op, Ap- 
pleton, WI, was used for all standard experi- 
mental test solutions. Rennet extract and 
Marzyme (single strength) containing 100 
rennet units/ml (RU/ml) were commercial 
milk-clotting preparations obtained from Miles 
Laboratories. Porcine pepsin (20,000:1) was 
purchased from Calbiochem as a dry powder. It 
was dissolved in .1 M KH 2 P04 to make solu- 
tions for assays. Other chemicals were reagent 
grade laboratory stock. 



2 100 - 




PRIMARY REACTION TIME (t, ((seconds) 



Figure 2. Graphical representation of equation 
[18]. 
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The assay substrate solution for all experi- 
ments was prepared in the following manner 
unless otherwise noted. For a final substrate 
volume of 100 ml, .4 ml of 1 N NaOH was 
added to 50 ml of .1 M KH2PO4 solution in a 
150-ml beaker. Five ml of .1 M CaCXi and 35 
ml of distilled, deionized water were mixed in 
another beaker. The two solutions were com- 
bined, and 2.00 g of dry skim milk was added 
to the combined solution while it was agitated 
with a magnetic stirrer. After the dry powder 
was dissolved completely, the solution was 
transferred to a 100-ml volumetric flask, 
diluted to the mark with washings from the 
beaker, and mixed thoroughly by repeated 
inversions. The substrate solution then was 
allowed to equilibrate for at least 2 h in the 
30 C water bath where clotting times were 
determined. 

Clotting times were determined in basically 
the same way as in the Berridge method (2) 
using a Sommer-Matsen rennet tester device. To 
accommodate 20-ml test tubes instead of the 
rennet bottles usually employed with the 
apparatus, the device was modified by #4 
cylinder guards added to the rollers. During 
clotting tests the tubes were rotated at 45 rpm. 
To obtain more consistent clotting times with 
smaller random error, five 4-mm glass beads 
were added to each assay tube. Assay tubes 
containing glass beads were preincubated in the 
water bath before use. 

For clotting time determination, 1 ml of a 
milk solution was pipetted into a prewarmed 
assay tube with a warm pipette. The tube was 
placed on the rollers of the rennet tester so that 
the lower half of the tube was submerged under 
the bath water, and the rotation was started. A 
predetermined amount of standard enzyme, 
usually 2 x 10~ 2 RU per sample, was injected 
directly into the assay milk just as a stop watch 
was started. The tube was shaken gently for 3 
to 4 s to disperse the enzyme; then it was 
replaced on the rollers. The watch was stopped 
at the first sign of coagulation in the tube. The 
elapsed time for each sample was recorded as 
its "secondary clotting time". 

Enzyme concentrations in an unknown 
sample were determined as follows. A known 
amount of an unknown sample was injected into 
a prewarmed assay solution with a 100-/il sy- 
ringe, and a primary timer was started to ascer- 
tain elapsed time from addition of enzyme. The 



solution was agitated thoroughly on a vortex 
mixer for 15 to 30 s, then replaced in the water 
bath. One-milliliter samples of the assay milk, 
now containing some enzyme, were transferred 
to assay tubes containing glass beads as de- 
scribed. At intervals, these solutions were used 
to make secondary clotting time determinations. 
The time between mixing of the assay solution 
with the unknown sample and the start of the 
secondary clotting time determination is 
referred to as the "primary reaction time." The 
relationship between the primary reaction time 
and secondary clotting time was used to deter- 
mine the enzyme concentration in the assay 
solution. This in turn could be related to the 
enzyme concentration in the unknown solution. 
A typical sample calculation is described in the 
Appendix. 

RESULTS 

Experimental Verification of Equation [18] 

The relationship between primary reaction 
time (ti) and secondary clotting time (t a ) was 
tested experimentally for several known primary 
enzyme concentrations with rennet extract. 
Figure 3 shows a comparison of clotting data 
and theoretical predictions of equation [18) 
for three concentrations (E0 of enzyme. There 
were differences between theory and data that 
cannot be attributed to experimental error. 
Specifically, this figure indicates that the 
secondary clotting times obtained at long 
primary reaction times are greater than those 
predicted by equation [18]. If these data are 




" 600 1200 1800 2400 3000 3600 4200 4800 5400 6000 
PRIMARY REACTION TIME (t,l (seconds) 

Figure 3. Theoretical and experimental relation- 
ship between the primary reaction time and secondary 
clotting time. (Assay solution clotted with rennet 
extract. The arrow indicates useful data range for 
accurate regressions for E,/E 2 = .010.) 
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used to calculate primary enzyme concentrations 
via slopes, low concentrations are obtained 
invariably. This problem also was observed for 
other clotting enzymes in this study. 

Careful examination of a large number of 
data showed that in the limit as ti the 
theory and experimental data coincide. Further- 
more, for all practical purposes, theory and 
experiment are equivalent for a significant 
range of t t . If assays are limited to the early 
portion of the t! vs. t a relationship, slopes that 
allow accurate determination of the primary 
enzyme concentration (Ej) can be obtained. 

Assays of solutions containing known 
enzyme concentrations (Ei) were made for 
rennet extract, pepsin, and "Marzyme". In all 
cases 2 x 10~ 2 RU/ml was the concentration 
of the secondary clotting enzyme (E 2 ). Con- 
centrations of Ej were varied to give dilu- 
tions (r) of .100, .010, and .001. The results 
from these assays are in Table 1 . The systematic 
error apparent at the highest primary enzyme 
concentration (r=.100) is due to inability of 
data to be taken at sufficiently small tj. Figure 
4 shows a more complete set of data for rennet 
extract. 

Use of Technique for Determination 
of Enzyme in Milk and Whey 

A sample of reconstituted skim milk (10% 
solids, .01 M added CaCl 2 , pH 6.6) was co- 
agulated by adding 2 X 10 — 2 RU of rennet 



extract per milliliter of milk. After enzyme 
addition, but before coagulation, 100 jul of this 
milk was added to 9.9 ml of assay solution. 
The enzyme concentration was determined by 
the methods described. The amount of enzyme 
from this assay technique was 2.04 ± .14 
X 10 -2 RU/ml (see Table 2). 

Fifty minutes after rennet addition, the milk 
gel that formed was cut repeatedly with a 
spatula, and the resulting slurry was centri- 
fuged at 900 X g for 10 min to separate whey 
and curd. The whey was analyzed for enzyme 
as described and contained 1.75 ± .08 X 10 — 2 
RU/ml. For the distribution of whey and curd 
(37.5 g whey, 12.5 g curd), this corresponds to 
66% recovery of the enzyme in the whey 
fraction, in agreement with observations of 
Holmes et al. (15). 

An amount of enzyme corresponding to a 
concentration of 2.00 X 10~ 2 RU/ml was 
added to a portion of this whey, and this 
sample was analyzed for enzyme content. The 
new concentration of enzyme, including the 
added enzyme, was 3.86 ± .14 x 10~ 2 RU/ml. 
Thus, the amount added was determined 
experimentally as 2.11 ± .16 X 10 -2 RU/ml, 
indicating that the method can determine 
enzyme accurately in whey solutions. Table 2 
summarizes these findings. 

Assay of Enzyme in Curd 

An attempt was made to assay directly for 
enzyme contained in curd by adding frozen and 



TABLE 1 . Determination of various enzymes. 



Calculated Slope 
Dilution factor 1 slope standard 

Enzyme r b deviation % Recovery 



Rennet extract 


.1000 


-.08179 


Porcine pepsin 


.1000 


-.08429 


Marzyme 


.1000 


-.07859 


Rennet extract 


.0100 


-.01025 


Porcine pepsin 


.0100 


-.00956 


Marzyme 


.0100 


-.00974 


Rennet extract 


.0010 


-.00098 


Porcine pepsin 


.0010 


-.00106 


Marzyme 


.0010 


-.00104 



SD 



.00115 


89.1 


1.2% 


3 


8 


.00110 


92.1 


1.1% 


3 


8 


.00082 


85.3 


.8% 


3 


8 


.00102 


103.5 


10% 





12 


.00118 


96.5 


12% 





12 


.00120 


98.4 


12% 





12 


.00007 


98 


7% 


8 


24 


.00011 


106 


11% 


8 


24 


.00008 


104 


8% 


8 


24 



'Dilution factor = concentration of primary enzyme/concentration of secondary enzyme = E,/E 2 where 
E 2 = 2 X 10 — 2 rennet units/ml in all runs, n+ = number of intercept points used in the determination, and 
m+ = number of nonintercept points used in the determination. 
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Figure 4. Use of the assay technique to determine 
the enzyme content of standard solutions containing 
rennet extract. 



ground curd to an assay solution. This technique 
yielded an assay of only 1.88 X 10~~ 2 RU/g of 
curd compared to an expected 2.75 X 10 — 2 
RU/g. This result is probably due to an in- 
complete release of the enzyme from the curd 
matrix. 

Other Assay Solutions 

The assay technique described should be 
expected to work with many assay solutions. A 
solution containing 2% wt/vol dry skim milk 
solids, .0125 M CaCl 2 , and .04MHC1 was used 



for some determinations on known samples. 
This assay solution produced satisfactory 
results for rennet extract and Marzyme deter- 
minations but slightly low results for pepsin 
determinations. This phenomenon was not 
investigated thoroughly or resolved, but there is 
some indication that the low pepsin deter- 
minations may have been due to nonspecific 
adsorption of the pepsin on the glass in the 
assay system. 

DISCUSSION OF RESULTS 

The experimentally observed result that the 
relationship between primary reaction time and 
secondary clotting time is not linear over the 
entire primary reaction range is unfortunate 
because it limits precision of this technique and 
makes it less convenient to use. Computer 
simulation studies have indicated that the 
deviation from theory seen in the experimental 
results could be caused by a nonenzymatic 
reaction, which must occur before clotting can 
be observed. Such a reaction likely involves 
hydrolysis products of K -casein. Because the 
physical flocculation of micelles involves this 
nonenzymatic reaction, experimental observa- 
tions can be explained partially in these terms. 

However, even with the inconvenience 
caused by this nonlinear relationship, this 
technique still is used more easily and more 
rapidly than other techniques. Enzyme con- 
centrations as small as 2 X 10 — 5 RU/ml in 
assay solutions are determined easily to ± 
10% precision in less than 2 h. This can be 
compared with a conventional assay limit of 
10 -4 RU/ml (which takes up to 48 h) in 



TABLE 2. Determination of enzyme in milk and whey. 





Amount added 


Amount determined 




Sample 


X 10~ 2 (RU'/ml) 


X 10~ 2 


(RU/ml) 


% Recovery 








SD 


SD 


Milk 


2.00 


2.04 


.14 


102 7 


Whey 




1.75 


.08 


88 4 


Whey + 2.00 X 10~~ 2 










RU/ml 


3.75 ± .08 


3.86 


.14 


103 4 


Difference in 










whey samples 


2.00 


2.11 


.16 


105.5 8 



1 Rennet units. 

2 Concentration of enzyme in whey relative to that in milk. 
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casein agar diffusion tests. Even considering 
that the unknown usually is diluted into the 
assay milk, it should be possible to assay 
solutions containing 10 — 4 RU in a much 
shorter time. The technique should be partic- 
ularly useful in assaying for enzyme in milk 
solutions that have been passed through im- 
mobilized enzyme columns. A continuous 
decrease of secondary clotting time after the 
solution exits the column would be an indication 
that there is soluble enzyme in the solution. 
The rate of decrease could be used to estimate 
the amount of soluble enzyme. 

It should be possible to use this general 
procedure in assaying for other enzymes with 
other assay solutions. In any enzymatic reaction 
system, the basic mathematical technique 
described herein can be used. For example, the 
enzymatic hydrolysis of benzoyl arginine 
ethylester by trypsin results in the release of a 
hydrogen ion. From a stock solution, the total 
conversion can be indicated by the change of 
pH that results from the reaction. By adding an 
indicator such asphenophthalein to the reaction 
mixture, a given extent of reaction can be 
indicated. The color change reaction in this case 
would be analogous to the clotting reaction. 
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APPENDIX 

Typical Example 

Three replicates of an assay solution con- 
taining 2 X 10~~ 5 RU were prepared by 5 fJ.1 of 
standard chymosin (.4 RU/ml) injected into 
each of three test tubes containing 10 ml of 
assay milk. One tube containing assay solution 
but no added enzyme was used as a control and 
for intercept determinations (see following 
discussion). All tubes were preincubated in a 
30°C water bath. 

After the 10-ml tubes were agitated thor- 
oughly, 1 ml samples of the solution were 
transferred with a warm pipet into prewarmed 
test tubes containing glass beads. At 10-min 
intervals an additional 2 X 10 — 2 RU/ml was 
added to these 1-ml subsamples by 5 /xl of 4.0 
RU/ml enzyme standard solution injected into 
the test tubes while they were rotated in the 



bath. Clotting times were determined with a 
stopwatch by the exact time measured from the 
secondary injection to the onset of visible 
flocculation. The experimental results for a 
control and three replicates are in Table 1-A. 

A significant improvement of precision of 
the determinations was obtained by inclusion in 
the regression analysis of experimental data from 
the secondary clotting times of control samples. 
Strictly speaking, these data do not belong 
to the regression set because for these deter- 
minations the total concentration of enzyme 
is only (E 2 ) rather than (Ei + E 2 ) as for the 
regular regression points. 

The secondary clotting time measures from 
control samples however, can be adjusted to 
represent intercepts (t! = 0) for unknown 
determinations by suitable modification of the 
experimental measures. As was pointed out, the 
intercept of the t] vs. t a plot equals the clotting 



TABLE 1-A. Experiment data for the example. 



Secondary clotting time (t a ) (s) 
(min) Control 1 



Primary reaction _ .. ^ 

. ' . Replicates 
time (t, ) ! 



10 


75.65 


74.94 


74.05 


74.31 


20 


76.36 


72.76 


73.36 


72.50 


30 


75.25 


72.99 


73.08 


72.41 


40 


75.72 


71.87 


72.97 


72.64 


50 


74.44 


72.55 


71.34 


71.59 


60 


74.11 


71.39 


72.48 


71.11 


70 


74.41 


70.42 


70.52 


69.97 


80 


73.81 


69.33 


69.93 


71.10 



'The data from replicates 1, 2, and 3 were pooled and fit to equation [18] . They were analyzed in terms 
of standard linear regression formulae to determine values for the slope (b), intercept (t' c ), standard deviation 
(S X y), and standard deviation of the slope (Sfo) (33). Results are in Table 2-A. 
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TABLE 2-A. Regression analysis. 








Replicates only 


Replicates plus intercept 


Slope (b) 

Slope standard deviation S5 
Data standard deviation S„„ 
Intercept 

Number of regression pts. 


-.00091 
.00009 
.58 

74.51 

24 


-.00098 
.00007 
.67 

74.75 

32 



time of a control sample with a concentration 
of enzyme equal to (Ei + E 2 ). The clotting 
time of a control sample containing a con- 
centration of enzyme E2 is (t c ) and is related 
to that of a sample containing a concentration 
(Ei + E 2 ) by equation [20] . If to in equation 
[20] is set equal to zero, then equation [20] 
becomes: 



Ei +E 2 



When data are fit to equation [18], the 
estimated slope of the regression line (b) is 
given by: 



= -E l /(E L +E 2 ) 



so that: 



1 + b = E 2 /(E! + E 2 ) 
Equation (A-l) then can be written as; 



= t c (l + b) 



[A-2] 



[A-3] 



[A-4] 



For a set of t c j for the control sample, a 

corresponding set t' c j of intercepts for the assay 
sample can be calculated if b is known. The 
calculated b depends on the intercepts (tj 
= 0, t a = tq) used in the calculations. This fact 



suggests that an iterative procedure to find b is 
appropriate. 

A good first estimate of b, (b ), can be 
calculated from equation 18 without including 
any intercepts. This will allow calculation 
of a first set of intercept values, t c (l)j according 



'(1) 



j - t CJ (l + b ) 



[A-5] 



[A-l] i m p rovec i estimate of b, b(\) then can be 

calculated by t' c (^h included in a second 



'1 

regression calculation by equation [18]. 
A revised estimate of t c , t c (2) can be made by 
using the new estimate of b, b(j) according to: 

^c' (2) - } = x °cP +b (l)> I A " 6 1 

and so on. 

This iterative procedure will converge to a 
constant b, usually in 2 to 3 iterations. This 
provides the most precise estimate of b that 
can be obtained with the data. 

Regression of the data (including the eight 
intercepts in Table 1-A) used this procedure. 
The regression results for such a calculation 
are in Table 2-A. 

A simple calculation allows one to determine 
the enzyme concentration (Ej) from slope 
data. In the case described here Ei = 1.96 ± .14 
X 10~ 5 RU/ml as opposed to the known 2.00 
X 10~ s RU/ml. 
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